Classical Information-Theoretical View of Physical Measurements and Generalized 

Uncertainty Relations 
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General characterizations of physical measurements are discussed within the framework of the 
classical information theory. The uncertainty relation for simultaneous measurements of two physical 
observables is defined in this framework for generalized dynamic systems governed by a Sturm- 
Liouville type of equation of motion. In the first step, the reduction of Kennard-Robertson type 
uncertainties due to boundary conditions with a mean-square error is discussed quantitatively with 
reference to the information entropy. Several concrete examples of generalized uncertainty relations 
are given. Then, by considering disturbance effects, a universally valid uncertainty relation is 
investigated for the generalized equation of motion with a certain boundary condition. Necessary 
conditions for violation (reduction) of the Heisenberg-type uncertainty relation are discussed in 
detail. The reduction of the generalized uncertainty relation due to the boundary condition is then 
tested experimentally by re-analyzing data for measured electron densities in a hydrogen molecule 
encapsulated in a fullerene Ceo cage. 



INTRODUCTION 

The problem of observation in quantum theory is a 
long-standing one, raised immediately after quantum me- 
chanics was first established. One of the most important 
topics is interpretation of the uncertainty relation, which 
was first expounded by Heisenberg[lj. He introduced his 
uncertainty relation as a principle of quantum mechan- 
ics through a Gedanken experiment regarding position 
and momentum measurements of a point particle using 
an imaginary gamma-ray microscope. KennardQ and 
Robertson^ generalized Heisenberg's uncertainty rela- 
tion and proved it mathematically as a relation between 
standard deviations (square roots of variances) associ- 
ated with a pair of non-commuting operators. Following 
the publication of their proof, there remained some con- 
fusion about the variance of distributions, mean-square 
measurement errors, and disturbances because of mea- 
surements. Recently, uncertainty relations have been 
clearly understood [J, [5|, and a universally valid uncer- 
tainty relation has been obtained 0-0] . However, it is 
still necessary to pursue a better understanding of un- 
certainty relations and their applications because the ob- 
servation problem is a key issue of quantum mechanics 
that has not yet been completely solved. Experiments 
claiming to demonstrate a violation of the Heisenberg- 
type uncertainty relation^ must be carefully considered 
from various points of viewfioj. 

In this study, we would like to propose a universal view 
of physical measurements based on the classical informa- 
tion theory. Here "classical" means that aspects such 
as commutation relations, Hilbert spaces, and operators 
on those spaces are not introduced a priori, as in con- 
trast to the usual discourse of the quantum measurement 
theory llj, [l2| . Hence, all concepts related to physical 



measurements are defined in terms of classical informa- 
tion theory. Probabilistic aspects of measurements are 



introduced through random variables in the framework 
of the probability theory. Solutions of a classical equa- 
tion of motion, such as the equation for charged spinor 
fields, represent a classical charge distribution. The rea- 
son why we restrict to the classical theory is that some 
counterpart of the "quantum effect", i.e., uncertainty re- 
lations, already arises in the classical field theory without 
any explicit quantization. One aim of this investigation 
is to clarify the boundary between classical and quantum 
mechanics in order to achieve a better understanding of 
the latter. 



PHYSICAL MEASUREMENT 

We define terms associated with a physical measure- 
ment according to the classical information theory (l3j| as 
follows. Let if be a random variable for a given physi- 
cal system described by the Y-tuple, 8 = {9i, ■ ■ ■ ,8n}, 
where Oi is the i th physical parameter. The set of all 
possible values of (9; e K, denoted by O, is called the 
parameter set. The random variable 9£ is distributed 
according to the probability density function i(x; 6) > 0, 
which is normalized as J^en ^ x K x 'i @) = ^ where x G R 
is one possible value of the whole event (= fi) . For phys- 
ical applications, we introduce the probability amplitude 
defined by 

|#r;0)| a =f(a:;0). 

It is assumed that the behavior of the physical system is 
determined by the probability amplitude rather than the 
probability density. This is by analogy with the quantum 
mechanical amplitude. However, the introduction of an 
amplitude does not immediately mean that the theory 
has been quantized. Considering the probability density, 
the probability amplitude is only by multiplication with a 
unitary operator, which implies that the amplitude con- 
tains more information than the density. 
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A peace of experimental apparatus is assumed to output 
numbers distributed according to the probability den- 
sity. Any resulting set of numbers, 5£ n — {x±,- ■ ■ , x n }, 
drawn independently and identically distributed (i.i.d.), 
is called the experimental data. The estimate the physi- 
cal parameter is called a measurement. Because the ex- 
perimental data is i.i.d., the corresponding probability 
density function can be expressed as a product: 

n 

f(ar n ;0) = JJ Kxj-,6). 

3 = 1 

A function mapping the experimental data to one possi- 
ble value of the parameter set such as 

Ti : JT„ -s- 8 : {x 1} • • • , x n } i-> 6 t 

is called an estimator for the i th physical parameter, 
denoted by Ti ( 3£ n ) = 0j . The experimental error in the i 
th physical parameter is defined as the root mean square 
error: 

ti = E[(Ti(%~ n ) - Oif] 1 / 2 , 

where 6i is the true value of the i th physical parame- 
ter. The true values of physical parameters are unknown, 
but a mean-square error can be reduced below any de- 
sired value by accumulating a sufficiently large amount 
of experimental data, thanks to the law of large num- 
bers. If the mean value of the experimental error is zero, 
Ee t [Oi — 9i] = 0, after accumulating an infinite number of 
statistics, that estimator is called an unbiased estimator. 
Among such estimators, the one giving the least error is 
called the best estimator. 

Simultaneous measurements of two physical parameters 
are described as follows. For the random variable, 3£ , 
with two physical parameters = {6i,6 2 } T , the ex- 
perimental data with n samplings can be expressed as 
2t n = {si,-- - ,x n } = {{x\x 2 )1,--- ,{x\x 2 )l}. The 
probability density of the experimental data is expected 
to be a Gaussian distribution with two variables: 



■exp(-(xj -ff) T a~ x (xj -0)/2, 



2tt|<t| 1 /2 

thanks to the central limit theorem. Here, er is the co- 
variance matrix of the data. When two parameters can 
be measured independently, the probability density func- 
tion of the data becomes 



f( Xj ;e)= [] 



1 



■ exp 



jx) - Ojf 
2aJ 



In this case, it is known that the best estimators for phys- 
ical parameters and experimental errors are given as fol- 
lows: 



T-j ( 3 ~n ) 



1 ™ 
n ^— ' 

1 n 

n ' 
i=i 



(1) 



for i=l,2. 

We now introduce quantitative informational proper- 
ties. The TV-tuple of random variables, 

<91ogf(x,0) 



Vi 



36,. 



is called the score. It can be shown that the expectation 
value of the score is zero [3]. The score represents the 
sensitivity of the experimental data to the i th physical 
parameter. A large value of the score means that the ex- 
perimental data is sensitive to the i th parameter and is 
expected to give a small error. The score for the experi- 
mental data can be shown to be, Vi{3£ n ) — Sj=i 
because the data is i.i.d. The covariance matrix of the 



score 



J, 



„(0) = / dx f(>;0) 



dlogf(x;0) 01ogf(x;0) 



d6 l 



86, 



(2) 



is called the Fisher information matrix (FIM). FIM of 
the experimental data can be shown to be Jij(3^ n ; 0) = 
nJij(9), once again, because the data is i.i.d. T = 
{Ti, • • • ,Tjy} is an unbiased estimator, and S(0) is the 
covariance matrix of the data. Then, the Cramer-Rao 
inequality 

no)>j-Ho) (3) 

holds [Hj true as a matrix inequality, i.e., each element 
on the left-hand side is greater than or equal to each 
corresponding element on the right-hand side. This is one 
of the ke y th ings to consider when discussing uncertainty 
relations 16 1. 



EQUATION OF MOTION 

The equation of motion (EoM) is a differential equation 
describing the time evolution of the physical system. We 
assume that the physical system is governed by the prob- 
ability amplitude, which is a solution of EoM. Here, we 
assume that EoM satisfies the following conditions: 1) It 
is a separable equation with respect to the time vari- 
able: 2) The spatial part after time-variable separa- 
tion will be a holomorphic Sturm-Liouville function (l7|: 
and 3) The Sturm-Liouville differential operator (SL op- 
erator) is self-adjoint. The SL-operator is defined by 

mx)] = A. Ux)^{x) 



q(x)ip(x) + Xr(x)tp(x), 



where p(x), q(x), and r(x) are smooth functions given by 
the dynamical system under consideration. Solving EoM 
leads to an eigenvalue problem of the form L[tf>(x)] = 
with the following boundary conditions on a finite inter- 
val [a, b]: 

dip(x) 



ip(a 
4>{b) + k 2 



dx 
dipjx) 
dx 



0. 
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Self-adjointness of the SL operator ensures that its eigen- 
functions form a complete orthogonal system. By using 
the complete orthogonal system obtained from EoM, the 
generalized Fourier transformation (GFT) can be defined 
by 



/(a)= / duiomv, 



(4) 



where £) is the eigen-function of EoM with eigen- 
value A. The existence of GFT with an appropriate inte- 
gration measure and their inverse transformations are en- 
sured by the generalized expansion theorem[IlS|. For 
GFT, the conservation of normalization of the two func- 
tions is ensured by the Parseval theorem 22j, |23j . This 
set of observables {£, A} is called the GFT dual pair. In 
the next section, we will obtain a non-trivial restriction 
on the uncertainty of a simultaneous measurement of the 
GFT dual pair. 



GENERALIZED UNCERTAINTY RELATIONS 
WITHOUT DISTURBANCES 

Equation of wave motion 

The equation of wave motion consists s typical example 
of EoM discussed in the last section. After separating 
out the time component, the steady-state solution in an 
infinitely large box is a plane wave solution of the form 

= aiexp + "2 exp ( 5 ) 

where £j are two physical variables (observables) with 
appropriate dimensions and h is a dimensional physical 
constant canceling out the dimensions of ^1^2 • It is well 
known that Eq. specifies a complete orthogonal sys- 
tem, and the integral transformation with a kernel as 
given by Eq. ([5]) is the usual Fourier (inverse Fourier) 
transformation. Suppose that the experimental appara- 
tus has been prepared such that the initial probability 
density for physical observables (£1,2) is the Gaussian dis- 
tribution 



f(&j {Mij 0*}) 



f 



2na, 



exp 



2a? 



(6) 



as would naturally be expected from the central limit 
theorem. Then, we assume that these two observables 
constitute the GFT dual pair. In the wave function 
case, GFT is equivalent to a usual Fourier transforma- 
tion. Probability amplitudes describing the above prob- 
ability densities are introduced as follows: 




x exp 



mi) 2 . .M6 - m) 
— h % 



The Gaussian distribution of Eq. © is obtained by 
squaring it. After performing GFT, the probability am- 
plitude of the GFT dual pair, £ 2 , becomes 




d£i V>i(£i;{At,""})exp 



h 



(6 - 

h 2 



1 — - — 



(7) 



Because the transformation kernel in Eq. ([5]) is a gen- 
eral solution of EoM as a function of the variable £1 , the 
GFT dual pair amplitude, ■02, is the solution in the GFT 
dual space, by direct analogy with the relation between 
configuration and momentum spaces. By squaring the 
probability amplitude of Eq. ([7]), one obtains the Gaus- 
sian distribution in the form 



f(6;{M2,o-2 = — }) 



2 1 

— =- exp 

7T h 



2(6 - V2?°\ 

h 2 



which has a standard deviation of a-i = h/2oi. Then, 
standard deviations of the two GFT dual parameters sat- 
isfy the relation 



0102 = 2- 



(8) 



Next we show that this relation gives a lower bound on 
measurement errors. Preparation of the physical system 
with observable £1 and measurement of the observable 
£2 for the same system can be performed independently. 
This is justified because we are treating a classical field 
that has the wave equation as EoM. Then, that total 
probability density can be expressed as f(£; {fi, <x}) = 
f(£ii {Mij a i})K^2i {^2, 02})- FIM can then be obtained 
from the above formula as follows: 



J{^~n] A*) 




d£,2 f(£;/i 



nJi(n). 



<91ogf(£; n) <91ogf(£;/x) 



(9) 



Here we omit the cr-component of FIM, because it is not 
related to the uncertainty relation. On the other hand, 
the covariance matrix for simultaneous measurements of 
the two independent parameters is given as follows: 



E(M) = 




(10) 



according to Eq. ([T]). By comparing Eqs. ([9]) and (|T0|) . it 
can be observed that this experiment gives the lower limit 
of the Cramer-Rao inequality [Eq. ©]. Then, Eq. (0 
gives a lower bound on the initial distributions of the two 
parameters, and this in turn yields Kennard-Robertson 
type uncertainty relation, a±a2 > h/2. 



4 



Hydrogen atom 

As a second example, we consider the following EoM 
of the SL type operator [24). 



dr 2 



rR(r) H rR{r) - K 2 rR(r) = 0. (11) 



This is the radial component of the Schrodinger equation 
in a Coulomb potential when the angular momentum is 
zero. If we consider a hydrogen atom, as — h 2 /m e e 2 — 
0.592 x 10 _8 cm is called the Bohr radius. It specifics a 
typical atomic length scale. Because we do not require 
any quantization, this is the equation for a classical elec- 
tron field called the de Broglie field. Solutions of this 
equation give the charge density distribution of a classi- 
cal electron field. A set of solutions normalized to unity 
in < r < 00 is given by 



Rn(r) 



( 2 V 1 



yasn J 2n 2 



■ exp 



' \ r i / 2r 



dsn 



dsn 



where L\(») are Laguerre polynomials. This set consists 
a complete orthonormal system satisfying 



dr r 2 R n (r)R m (r) = S n 



(12) 



The corresponding eigenvalues can be obtained as E n — 
— h 2 K 2 /2m e with k — (oBfi) -1 . Once again we take the 
initial probability density of the electron to be a Gaussian 
distribution: 



ip(r;a r ) = 



*k) ' (13) 



f(r;ay) = \ip(r; a r )\ 2 , 



( 1 



V27ra r 2 ; 6XP V 2al 



(14) 



One reason why we start with the Gaussian distribution 
is that it is expected to be obtained from the central limit 
theorem for the physical measurements of i.i.d. and will 
lead the minimum uncertainty condition. GFTs of these 
eigen-functions are given as 

/>00 

<j> n {<r r ) = C(a r ) / drr 2 i>(r;<j r )R n (r), (15) 
Jo 

= C (a r ) J] kz l "> m2l+m ( m + 2 ) r W 



x F 



^„ 7T4 m\T{n — m)(aon) 2 



m + 3 1 / ay 



2 2 \aon 




(16) 



where _F(»,»;») is the confluent hypergeometric series, 
and C(ay) is an appropriate normalization factor. Be- 
cause the solutions Ri{r) form a complete orthonormal 
set of functions, as shown in Eq. (|12"j) . the Gaussian 



distribution can be expressed as the inverse GFT of an 
infinite summation: 



00 

V>(r;oy) = }^4>i(a r )Ri(r). 

z=l 



We have checked numerically that this inverse GFT 
would transfer <^>j(ay) back to the original Gaussian dis- 
tribution, except that it would be close to the origin. 
Note that atomic units (a.u.)[25| are used in all numer- 
ical calculations in this study (m e = 1, h — 1, and 
e = 1). In this case, the usual Fourier transformation 
of ip(r; ay) does not have any clear physical meaning be- 
cause the plane wave solution has no vanishing values 
at r = 00, whereas tj)(r — > 00; ay) = for a hydrogen 
atom. The GFT dual parameter of the radial coordinate 
r must be an energy eigenvalue E n . Probability density 
functions in the GFT dual space have Gaussian distribu- 
tions, as shown in FIG. [T] The relation between stan- 
dard deviations of the initial Gaussian distribution and 
the corresponding energy spectrum is shown in Fig. [2] 
We can observe that the uncertainty of the energy spec- 
trum decreases position uncertainty increases. Moreover, 
products of these two uncertainties are almost constant 
at <r r <TE — 1/1-72. This condition also gives the lower 
bound from Cramer-Rao inequality because the GFT 
dual probability density function (<p n ) has a Gaussian 
distribution, and the discussion from the first example 
can still be applied. A generalized uncertainty relation 
for this case can be expressed as 

° r ' aE - T72"' ( > 

where E a is the atomic unit of energy. This relation 
is similar to the Kennard-Robertson uncertainty, except 
that the coefficient in front of the dimensional parameters 
on the right-hand side is 1/1.72, instead of the usual value 
of 1/2. 

Let us look at this new relation from another point of 
view. The reason why a hydrogen atom does not collapse 
to the size of the proton is usually discussed with only 
a qualitative reference to the uncertainty relation. Here, 
we investigate this matter in a quantitative way using 
the new uncertainty relation. The total energy of the 
electron can be given as 

p 2 1 

Eh — 3— , 

2 r 

in a.u. The first and second terms are basically ki- 
netic and potential energies. The factor of 3 in front 
of the kinetic energy originates from the spatial degree 
of freedom. Suppose now that the uncertainty relation is 
<J p (T r > f u . In a ground state hydrogen atom, the elec- 
tron is confined within the Bohr radius, ay = «b = To, 
which suggests that oy = f u /cr r — f u /ro- The total en- 
ergy can thus be written as 

f 2 1 

J u_ 



n 



a 

r 



FIG. 1. Example of the probability density function in the 
GFT dual space for the radial distribution of the classical 
electron field in a hydrogen atom. The solid line is a Gaussian 
fit. A Gaussian distribution is assumed for the radial variable, 
with variance = 10 2 . 



FIG. 3. Entropy obtained from hydrogen energy states and 
radial distributions as a function of the width of the radial dis- 
tribution (ay). The entropy of hydrogen energy states reaches 
a minimum at the Bohr radius. 
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FIG. 2. Relation between standard deviations of the ini- 
tial Gaussian distribution and the corresponding energy spec- 
trum. The electron mass and Bohr radius are set to unity 
(a.u.) in calculations. 
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FIG. 4. Width of the radial distribution of the beam mul- 
tiplied by its momentum as a function of the beam width. 
GFT can maintain the product of widths of the pair at 
(jROp = 1/0.30, where h is set to unity. 



where ro is the hydrogen radius in the ground state. The 
value ro must give the minimum energy, then, it can be 
obtained as a solution of dE^/dro = 0, i.e., r = 3/^. 
If we set ro = as = 1, then, we get j u = l/\/3 = 
1/1.732 • • • which is close to the factor in the new uncer- 
tainty relation. If one can assume, that the uncertainty 
relation gives the exact hydrogen radius, then the un- 
certainty relation between GFT dual parameters for the 
system under the SL operator of Eq. (jlip must be 



o> • cte > 



a-BE a 

~7T' 



(18) 



We can give yet another view of the hydrogen atom 
in terms of the classical information theory. The entropy 
of the energy spectrum of the hydrogen atom can be cal- 



culated numerically from Eq. (|16j) as 

oo 

S(E) = -} n (oy) log0„(oy). 

n=l 

Numerical results are shown in FIG. as a function of 
the width of the radial distribution. The global minimum 
point of the entropy coincides with the Bohr radius, as 
shown Fig. ([3]). 

Classical matter field in a cylinder 

We now consider a particle beam entering a cylindri- 
cal beam pipe. The Schrddinger equation in cylindrical 
coordinates is considered as EoM for the classical matter 
field. EoM for a radial variable R(r) with zero angular 
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momentum can be expressed as 



d 
dr 



dR(r) 
dr 



-R(r) + {k 2 - X 2 )rR(r) = 0, 



which is also an SL type equation. Solutions with the 
infinite potential barrier at r = r are obtained as 



Rf(r) 



V2 



ro\Ji+i(z?)\ 



J, 



ro 



(19) 



where Ji(x) is the I th order Bessel function of the first 
kind and zf is its n-th zero point. This series of functions 
constitutes an orthonormal system with normalization 
given as 



dr rR\\r)R\ n {r) = S n 



The corresponding eigenvalues are 



fc 2 = — 
\ r o 

P? = hk, 

= h —. 



A 2 , 



(20) 



Here, A is the eigenvalue from the appropriate boundary 
condition in the ^-coordinate and will be neglected in the 
following discussion. Even though a physical constant, 
H, arises in the eigenvalues, this is a classical theory, and 
the solution i?"(r) represents a classical de Broglie field. 
Suppose a particle beam with Gaussian distribution in 
the radial direction of the cylindrical coordinate enter a 
cylinder of radius ro . The beam is assumed to be coaxial 
with respect to the cylinder, so the angular momentum 
of the beam is zero (1 = 0). This radial distribution can 
be expanded in terms of the solutions of Eq. (TT9l) as 



dr r ip(r;a R )RQ{r), 



where ip(r;a R ) is the radial distribution of the incident 
beam represented by Eq. (fT3|) . When a R <C r , the 
above integration can be performed analytically to give 



4ct 



3/2 



(27r)i/4 roJl ( z «) 



exp 



o-fl go 

ro 



The usual Fourier dual pair has no clear physical meaning 
because the solution of EoM with a cylindrical boundary 
condition is not a plane wave. The GFT dual pair can 
be given as {r, Pq = z^/r }, and GFT transforms the 
Gaussian to another Gaussian. The width of the Pq dis- 
tribution is proportional to the reciprocal of a R , as shown 
in FIG. H The relation between the widths of the GFT 
pair is found numerically to be (Jr&p = h/3.0, which 
gives a smaller uncertainty than the Kennard-Robertson 
relation expected from the usual Fourier dual pair. This 
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FIG. 5. Entropies of the radial distribution of the beam and 
those of the GFT dual pair of the momentum eigenstate dis- 
tribution, h and ro are set to unity. 



reduction of the uncertainty can be explained with re- 
spect to the entropy. The entropy of the radial distribu- 
tion, which is contained within the cylinder, can be cal- 
culated from the probability density function [Eq.(|14])] 

as 



b (r < r ; o-fl) 



dr f(r; a R ) logf(r; a R ), 



whereas the entropy of the momentum distribution can 
be obtained as 

oo 

S c y\p; a R ) = - £ <j> n {a R ) log0 n (a fl ). 

n=l 

As the beam width increases, the entropy of the radial 
distribution increases, whereas that of the momentum 
distribution decreases such that their sum is maintained 
to be constant, as shown in FIG. [5] This behavior of en- 
tropies is easily understood qualitatively if we recall that 
the entropy is a measure of the total amount of informa- 
tion in distributions. Under the minimum uncertainty 
condition, the total amount of information can be main- 
tained constant. Here, let us consider them quantita- 
tively in order to estimate the information gain obtained 
from the boundary condition. When there are no bound- 
ary conditions on radial distributions, the information 
entropy of the usual Fourier dual pair is 



-iGauss i 



1 



(a x ) + S^ ss (a p ) = - log (4ttV^) + 2A. 

where A is an arbitrary value due to ambiguity in the 
choice of integration measure. Then, the relation be- 
tween the entropy and mean-square errors can be ex- 
pressed as 



<t x <j v 



1 

27 

h 



exp(S(a x ) + S{a p )-l-2A), 



exp(S(a x ) + S(cr p )). 
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FIG. 6. Width of the radial distribution of the initial field 
multiplied its GFT dual momentum. GFT can maintain the 
product of widths of the pair at a T ap — 1/2.9, where h is set 
to unity. 



• Orthonormality relation: 



dr r 2 R n (r)R m (r) = 5„ 



Here, we set A equal to — 1/2 log (eirti) to obtain the min- 
imum uncertainty relation when S(a x ) + S(a p ) — 0. The 
average information due to the cylindrical boundary con- 
dition can be estimated as follows. The radial distribu- 
tion of the incident beam is assumed to be the Gaussian 
distribution, as given in Eq. ([6]). The total amount of 
information in the cylinder of radius ro can be calculated 
as 

/ro 
drf(r,{0,a R }) 
- r 




where erf(») is the error function. The entropy within 
the cylinder can be calculated as 

S cyl (r ■:<J R ) = -I(r 07 a R ) log I(r ,a R ). 

The maximum value of the total information is found to 
be S cyl = e _1 by solving dS cyl /dr$ — 0. This amount 
of information gain due to the cylindrical boundary con- 
dition can decrease the minimum uncertainty according 
to 




h 

~ 2J' 

which is consistent with the numerical result shown in 
FIG. |H This reduction of the uncertainty relation is due 
to the information gain from the cylindrical boundary 
condition. 



• Eigenvalue for the radial equation: 

P n = hk 

717T 

= K — , (n = ±l, ±2,---) 

ro 

• GFT and inverse GFT for the Gaussian distribu- 
tion: 

0n(<7R) = / dr r 2 ip(r;cr R )R n (r) 
Jo 

oo 

4>(r;<TR) =^2<t)i((TR)Ri(r) 

• GFT dual integration with approximation a R <C 

p OO 

<f>n{cr R ) ~ / dr r 2 ip(r;a R )Rn(r) 
Jo 

2 9 / 4 ^/ 4 5/2 ( rvKa R V 

• Generalized uncertainty relation(see FIG. [6]): 

h 



Classical matter field in a sphere 

The last example is a classical matter field confined in 
a sphere. The calculations are almost the same as those 
in last the section, we only list the results here: 



In this case, mean-square errors once again give the 
Kennard-Robertson type uncertainty relation with a fac- 
tor of h/2.9 instead of h/2, which is consistent with 
the information gain because of the boundary condition, 
same as that as in the previous example. 
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GENERALIZED UNCERTAINTY RELATIONS 
WITH DISTURBANCES 

The universally valid uncertainty relation 

Next we consider a disturbance due to the measure- 
ments themselves. It is assumed that the initial con- 
ditions of the physical system agree with the minimum 
uncertainty condition, with mean values of {/Ji,a 2 } and 
variances of {^ 2l <J 2 }- The two observables {ni,^ 2 } are 
assumed to be a GFT-pair. Estimators associated with 
this measurement are assumed to be unbiased. Af- 
ter simultaneous measurements of physical parameters 
{fix, [i 2 }, final distributions are expected to be Gaus- 
sian because of the central limit theorem. The distur- 
bance, 5f, is defined as an increase in variance in the 
distribution of the i th observable after the measure- 
ment. The convolution of two Gaussian distributions 
with mean values and variances as {/ii, a 2 } and {/j, 2 , of} 
yields the Gaussian distribution with {jUi, of = of + Sf} 
and {/i2, erf = o\ +<y|}- Initial and final distributions are 
not necessarily a GFT pair. To find the minimum uncer- 
tainty condition after the measurement, we assume that 
the measurement has been performed with very weak 
coupling and gives a minimum uncertainty pair of dis- 
turbances with (5i<52 = h/2. Then, standard deviations 
after the measurement can be written in the form 



^2-2 

a l a 2 



Min^of] 




(21) 



The minimum uncertainty condition after the measure- 
ment is given by Si as da\a\jd8\ — such that 

h 2 of 



= 25i 0% - 



=> 01 = 




= 0, 



Then, the formula 



0\0 2 > (Jl<T 2 + — 



(22) 



(23) 



follows from (1211) and ([221) . Under this condition, we de- 
rive the universally valid uncertainty relation (UVUR) @- 



C1S2 + #102 + 5\8 2 > 2' 



(24) 



From Eq. (|23p and the positivity of standard deviations, 
Cji^i > 0, one can obtain 

(01 + 6 1 f{a 2 + S 2 ) 2 > (of + 5 2 ){cr 2 2 + 5 2 2 ), 

~2~2 
= ^1^2, 



> Ci<7 2 



Then, the inequality 



o\o 2 + <7i<5 2 + 5i<7 2 + 5iS 2 > <7\cr 2 + — 

is obtained. UVUR, Eq. (|M| . follows immediately from 
this formula [26j. 



Example: Hydrogen atom in a Cgo cage 

We now discuss the possibility of observing a violation 
of the lower bound in the Heisenberg-type uncertainty 
relation, <5i<52 > ft/2 (h is written as h in this section.). 
We have re-investigated position measurements of elec- 
trons in a hydrogen molecule trapped in a fullerene Ceo 
cage. In 2005, Komatsu et al. succeeded in encapsulating 
molecular hydrogen in fullerene C™ (H 2 @Gqo) with an 
efficiency of approximately 100% 27]. Sawa et al. mea- 
sured the electron density in the closed Ceo cage using X- 
rays from a syncrotron light source[28j. The electron den- 
sity was measured using X-rays of wavelength 0.0998 nm. 
The H 2 @Cgq sample was maintained at 50K. It has been 
confirmed experimentally that the rotational mode of 
H 2 @Cqq is almost eliminated at this temperature[29| . 
The inner diameter of fullerene Cqo is known to be 
approximately 0.7 nm. Here, we employ a radius of 
0.36 nm, as used in [3(j. Electrons belonging to H 2 @Ceo 
can be described by the Schrodinger equation with the 
boundary condition that the electron density is zero for 
r > 0.36 nm. The equation for the hydrogen atom has 
been solved numerically using Mathematical^. Here, 
we assume that energy eigenvalues are well approximated 
by those of the hydrogen atom for smaller principal quan- 
tum numbers. The solution proportional to e +r , which is 
abandoned as a solution for molecular hydrogen in vac- 
uum, is allowed for H2 in the Cqo cage. Energy eigen- 
values are listed in TABLE. HI along with the values in 
vacuum. All solutions have zero angular momentum, be- 
cause there no solutions with non-zero angular momen- 
tum in a spherically symmetric cage of radius 0.36 nm. 
Before the measurement, the electrons must be in the 
ground state. The uncertainty in the electron energy 
originates only from thermal fluctuations in electron en- 
ergies, which is estimated tobecrE = 4.31xlO _3 eVata 



n 


H2(Vac.) 


H 2 <i 


3>C 60 (Sol.l) 


H 2 i 


SC 60 (Sol.2) 


1 


-13.62 




-13.62 




-13.62 


2 


-3.41 




-1.14 




-1.14 


3 


-1.51 








1.14 


4 


-0.85 








13.62 



TABLE I. Energy eigenvalues of a hydrogen atom with bound- 
ary condition from the C60 cage. Here, n is the principal 
quantum number. In the C60 cage, the solutions proportional 
to e +r (Sol. 2) are permitted as are those proportional to e~ r 
(Sol.l), in contrast with the situation in vacuum (Vac). 
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temperature of 50 K. The electron position in the ground 
state has a distribution with variance equal to the square 
of the Bohr radius. Hence, the Bohr radius is considered 
as the uncertainty in the initial electron position, i.e., 
oy = as = 5.29 x 10~ 2 nm. The electron may be ex- 
cited to one of the excited states listed in TABLE, fl] after 
the measurement. Then, the energy uncertainty of the 
electron is at the most Se = 27.24 eV. The uncertainty 
in the electron density measurement is not clearly given, 
but the distribution of the electron density with respect 
to the radius has been measured very clearly, and the 
estimated number in the electrons in H2@Gqq is given as 
1.9±0.4(a" 



32j. If we assume 5 r — 1 x 10~ 2 nm, the calcu- 
lated electron density distribution is consistent with the 
figure and the estimated number of electrons in the cage 
is found to be 2.0±0.4, which is consistent with the exper- 
imental result. We thus use the value S r = 1 x 10 -2 nm 



for the position uncertainty of the measurement 33]. In 
conclusion, we obtain the following uncertainties: 



Se x err 
<je x Sr 
Se x Sr 
at,E 



1.44 eV nm, 

4.31 x 10~ 5 eV nm, 

2.72 x 10 _1 eV nm, 

3.60 x 10 _1 eV nm. 



Here, we observe that Heisenberg uncertainty pair, 
5e x Sr, exhibits smaller values than the expected vales 
a b E /2. 



SUMMARY AND CONCLUSIONS 

We have formulated the physical measurement process 
on the basis of the classical information theory. All 
processes related to measurements are treated without 
introducing any quantization of the physical system. 
In other words, all physical quantities are described 
by "c- numbers" . The probabilistic behavior of physical 
measurements comes from the assumption that the 
experimental data are random variables obeying some 
probability low. Even in this classical context, essential 
properties of uncertainty relations are exhibited. The 
Kennard-Robertson type of uncertainty relation arises 
for two physical quantities related to one another 
by Fourier transformation. Although Fourier trans- 
formation plays an essential role, it is not the only 
transformation leading to an uncertainty relation. We 
have investigated GFT dual pairs of physical quanti- 
ties governed by the Sturm-Liouville type differential 
equations and have obtained a generalized uncertainty 
relation for such GFT dual pairs. In addition, we 
have shown that Gaussian distributions realize the 
minimum uncertainty condition using the Cramer-Rao 
inequality, and the minimum uncertainty condition can 
give a smaller error than that implied by the usual 



Kennard-Robertson lower limit. This reduction of 
uncertainties can be understood quantitatively in terms 
of the information entropy because of boundary condi- 
tions. Another important result of this study is that 
the probability amplitude, which is defined as a square 
root of the probability density, plays as essential role 
in introducing the GFT dual pair of physical quantities 
and obtaining the generalized uncertainty relation. 
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